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The physical properties of magnetic nanoparticles have been investigated 
with focus on the influence of dipolar interparticle interaction. For weakly 
coupled nanoparticles, thermodynamic perturbation theory is employed to 
derive analytical expressions for the linear equilibrium susceptibility, the 
zero-field specific heat and averages of the local dipolar fields. By intro- 
ducing the averages of the dipolar fields in an expression for the relaxation 
rate of a single particle, a nontrivial dependence of the superparamagnetic 
blocking on the damping coefficient is evidenced. This damping dependence 
is interpreted in terms of the nonaxially symmetric potential created by the 
transverse component of the dipolar field. 

Strongly interacting nanoparticle systems are investigated experimen- 
tally in terms of spin glass behavior. Disorder and frustration arise in sam- 
ples consisting of frozen ferrofluids from the randomness in particle position 
and anisotropy axis orientation. A strongly interacting FeC system is shown 
to exhibit critical dynamics characteristic of a spin glass phase transition. 
Aging, memory and rejuvenation phenomena similar to those of conven- 
tional spin glasses are observed, albeit with much weaker rejuvenation effects 
than in both a Ag(ll at% Mn) Heisenberg and an Feo.sMno.sTiOs Ising spin 
glass. Differences in the nonequilibrium dynamics of the strongly interacting 
nanoparticle system and the two spin glass samples are discussed in terms 
of anisotropy and different timescales, due to the much longer microscopic 
flip time of a magnetic moment than of an atomic spin. 
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1 Introduction 



Ferro- and ferrimagnetic nanoparticles are important examples of how a 
reduction in size changes the properties of a magnetic material. For small 
particles it is energetically favorable to avoid domain walls and form only one 
magnetic domain. The magnetism of such single-domain particles has been 
an active field of research since the pioneering work of Stoner and Wohlfarth 
[1] and Neel [2] in the late 1940s. Because of new fabrication methods and 
characterization techniques, understanding of and interest in nanosized ma- 
terials have increased explosively within the disciplines of physics, chemistry, 
material science, and medicine. This development is driven by a large num- 
ber of applications; nanosized magnetic materials are used in, for example, 
magnetic recording media, ferrofluids, catalysts, and refrigerators; as well as 
by a large interest of fundamental nature. Nanomagnets made up of a small 
number of spins can be used to study quantum tunneling of magnetization 
[3]. In ferrofluids the dipolar interparticle interaction can be tuned by the 
particle concentration, and frozen ferrofluids have been shown to change 
their magnetic behavior from superparamagnetic at low concentrations to 
spin- glass-like in dense systems. 

The research on spin glasses started in the 1970s after the discovery by 
Cannella and Mydosh [4] of a peak in the ac susceptibility of diluted gold- 
iron alloys. Several different materials with various interaction mechanisms 
were soon found to exhibit this "new" magnetic behavior, all with two prop- 
erties in common — disorder and frustration. Spin glasses have since been 
widely studied, partly because they are excellent model systems of mate- 
rials with quenched disorder. An understanding of spin glasses can thus 
contribute to the understanding of other, more complex disordered systems, 
such as ceramic superconductors, polymers, gels, and dense nanoparticle 
systems. 

This article reviews the dynamic properties of magnetic nanoparticle 
systems with different interparticle interaction strength. In section 2 we 
discuss basic properties of noninteracting particle systems, and thermody- 
namic perturbation theory is used to study weakly interacting particle sys- 
tems. In section 3 we discuss the behavior of strongly interacting magnetic 
nanoparticle systems in the light of recent results in the field of spin glasses. 

2 Single-Domain Magnetic Nanoparticles 

The study of single-domain magnetic particles has been an active field of 
research since the pioneering work of Stoner and Wohlfarth [1], who studied 
the hysteretic rotation of the magnetization over the magnetic-anisotropy 
energy barrier under the influence of an applied field, and Neel [2], who 
predicted that at nonzero temperature the magnetization can overcome the 
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energy barrier as a result of thermal agitation. Later, Brown [5] derived 
the Fokker-Planck equation for the probability distribution of spin orienta- 
tions, starting from the stochastic Landau-Lifshitz equation, and calculated 
approximate expression for the relaxation time of particles with uniaxial 
anisotropy. The theoretically most well studied systems are noninteracting 
classical spins (representing the magnetization of the nanoparticles) with 
axially symmetric magnetic anisotropy A great step forward in comparing 
experiments and theory was taken when measurements on individual parti- 
cles were reported [6] . A profound knowledge of the physical properties of 
isolated particles is a prerequisite for further studies of phenomena such as 
quantum tunneling in molecular nanomagnets or dipole-dipole interaction 
in dense samples. 

In this section, some general properties of magnetic nanoparticles are 
first recalled (for more details, see, e.g., Refs. [7-9]). The definition of a 
magnetic nanoparticle is rather wide and includes ferro- and ferrimagnetic 
materials (e.g., 7-Fe203, Fe304, and Fei- x C x ) as well as magnetic nan- 
oclusters (e.g., Mni2 and Fes). Subsequently, thermodynamic properties of 
spins weakly coupled by the dipolar interaction are calculated. Dipolar in- 
teraction is, due to its long range and reduced symmetry, difficult to treat 
analytically; most previous work on dipolar interaction is therefore numeri- 
cal [10-13]. Here thermodynamic perturbation theory will be used to treat 
weak dipolar interaction analytically. Finally, the dynamical properties of 
magnetic nanoparticles are reviewed with focus on how relaxation time and 
superparamegnetic blocking are affected by weak dipolar interaction. For 
notational simplicity, it will be assumed throughout this section that the 
parameters characterizing different nanoparticles are identical (e.g., volume 
and anisotropy). 

2.1 General Properties 

The current studies of magnetic single-domain nanoparticles are limited to 
systems where the particles are fixed in space (realized, e.g., in frozen fer- 
rofluids, single crystals of molecular magnets, and magnetic nanoparticles 
in a solid matrix). We will also assume that every single-domain nanoparti- 
cle is in internal thermodynamic equilibrium and that its constituent spins 
rotate coherently. Moreover, we are considering only temperatures much 
lower than the Curie temperature, so the spontaneous magnetization is ap- 
proximately constant with temperature. Hence, the only relevant degree of 
freedom is the orientation of the net magnetic moment. 

The Hamiltonian of a single isolated nanoparticle consists of the mag- 
netic anisotropy (which creates preferential directions of the magnetic mo- 
ment orientation) and the Zeeman energy (which is the interaction energy 
between the magnetic moment and an external field). In the ensembles, 
the nanoparticles are supposed to be well separated by a nonconductive 
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medium (i.e., a ferrofluid in which the particles are coated with a surfac- 
tant). The only relevant interparticle interaction mechanism is therefore the 
dipole-dipole interaction. 

2.1.1 Magnetic Anisotropy 

The term magnetic anisotropy is used to describe the dependence of the 
internal energy on the direction of the spontaneous magnetization of the 
ferro/ferrimagnetic nanoparticle, creating "easy" and "hard" directions of 
magnetization. In general, a bulk sample of a ferromagnet exhibits mag- 
netic anisotropy with the same symmetry as in the crystal structure. This 
anisotropy energy originates from spin-orbit coupling and is called magne- 
tocrystalline anisotropy [14]. The two most common symmetries are uniaxial 
and cubic. For uniaxial symmetry the energy is given by 

£™ ni = KiV sin 2 9 + K 2 V sin 4 9 + ■ ■ ■ , (1) 

where V is the particle volume, K\ and K 2 are anisotropy constants, and 
9 is the angle between the magnetic moment and the symmetry axis. For 
cubic symmetry the anisotropy can be expressed in terms of the direction 
cosines (oj) as 

^cubic = # lF ( a 2 a 2 + a 2 a 2 + a 2 a2) + ^ 2 y a 2 a 2 a 2 + . . . ^ (2) 

where the a» are defined through a.\ = sin 9 cos 4>, a 2 = sin 9 sin <p and Q3 = 
cos#, 9 is the angle between the magnetization and the z axis, and <p is the 
azimuthal angle. 

For a single-domain ferromagnet, any nonspherical particle shape gives 
rise to shape anisotropy due to the internal magnetostatic energy. The mag- 
netostatic energy, for an ellipsoid of revolution, is equal to 

E m = \ix Q VM*{N z cos 2 9 + N x sin 2 9), (3) 

where 9 is the angle between the magnetic moment and the polar axis z, M s 
is the saturation magnetization, N z is the demagnetization factor along the 
polar axis, and N x = N y is the demagnetization factor along an equatorial 
axis. Both the magnetostatic energy for an ellipsoid and the uniaxial mag- 
netocrystalline anisotropy energy [Eq. (1)] can, to first order — except for 
a constant term — be written as 

E a = -Acos 2 9, (4) 

where A = KV is the anisotropy energy barrier and the uniaxial anisotropy 
constant K = ^fioM^(N x —N z ) in the case of shape anisotropy. For a prolate 
ellipsoid, K > and the anisotropy is of "easy axis" type, since there exist 
two minima of the anisotropy energy along ±z (the anisotropy axis). For an 
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oblate ellipsoid, K < and the anisotropy energy has its minimum in the 
whole xy plane. In this case the anisotropy is of "easy plane" type. 

With decreasing particle size, the magnetic contributions from the sur- 
face will eventually become more important than those from the bulk of the 
particle, and hence surface anisotropy energy will dominate over the magne- 
tocrystalline anisotropy and magnetostatic energies. A uniaxial anisotropy 
energy proportional to the particle surface S 

E ^c = KsScQS 2 e (5) 

has been observed experimentally by ferromagnetic resonance [15]. 

Hereafter, we will assume uniaxial anisotropy, of easy-axis type, given by 
Eq. (4) (if not otherwise indicated), since it is the simplest symmetry that 
contains the basic elements (potential minima, barriers) responsible for the 
important role of magnetic anisotropy in superparamagnets. Experimental 
evidence for uniaxial anisotropy is given in Refs. [15, 16]. 

2.1.2 Superparamagnetic Relaxation 

The uniaxial anisotropy energy creates two potential wells separated by the 
energy barrier A. The magnetic moment is subjected to thermal fluctua- 
tions and may undergo a Brownian-type rotation surmounting the potential 
barriers. This relaxation process was proposed and studied by Neel in 1949 
[2] and further developed by Brown in 1963 [5]. In the high potential barrier 
range, (3 A 3> 1, where = l/k^T, the characteristic time for the overbarrier 
rotation r can approximately be written in the Arrhenius form 

r ~ r exp(/L4) , (6) 

where To ~ 10~ 9 — 1CP 12 s. For observation times i bs much longer than the 
relaxation time, m maintains the thermalequilibrium distribution of orien- 
tations as in a classical paramagnet; however, because of the much larger 
magnetic moment than a single spin, this phenomenon was called super- 
paramagnetism [17]. The condition of superparamagnetism (t Q b s 3> t) cor- 
responds to a temperature range that fulfills ln(i b s /To) > (3 A . For t ^ s ~ 10 
s, due to the small value of To, this equilibrium range extends down to low 
thermal energies compared to the anisotropy energy (25 > [3 A). Hence, 
within the equilibrium regime, the system displays an isotropic behavior at 
high temperatures ((3 A <C 1), but a strongly anisotropic behavior at low 
temperatures {(3 A 3> 1). 

If t bs <C t, the magnetic moment is blocked in one of the potential wells, 
a state that corresponds to stable magnetization in a bulk magnet. If the 
measurement time is of the same order as the relaxation time (t bs ~ t), 
dynamical time-dependent effects are observed. 
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2.1.3 Effects of a Magnetic Field 



The Hamiltonian of a noninteracting nanoparticle with uniaxial anisotropy 
is given by 

A 

H = ^(m ■ n) 2 — /j, rh ■ H (7) 



where m is the magnetic moment with m = M S V, A = KV and n is a unit 
vector along the symmetry axis of the anisotropy energy (anisotropy direc- 
tion). By introducing unit vectors for the magnetic moment (s = fh/m) and 
the external magnetic field (h = H/H) and defining dimensionless parame- 
ters for the anisotropy and magnetic field 

a = f3A, £ = (3[i mH, (8) 

we can write a dimensionless Hamiltonian as 

-l3H = a{s-nf + £(s-h). (9) 

The bistable character of the zero-field Hamiltonian will be destroyed by 
a sufficiently large field. The critical field for n\\h is called the anisotropy 
field, and is expressed by 

2A 2K 

Hk = = — rp 10 

Horn /jl M s 

We can define another dimensionless field quantity 

h = w K = l< < n) 

which is the field measured in units of the anisotropy field. The Hamil- 
tonian, as a function of the angle 9 between the anisotropy axis and the 
magnetic moment (s ■ ft = cos 6), is shown in Fig. 1 for different values of 
the longitudinal field. 



2.1.4 Interparticle Interaction 

Dipole-dipole interaction is present in all magnetic spin systems, but usually 
other interaction mechanisms such as exchange interaction dominate. The 
relative weakness of the dipolar coupling between magnetic ions in param- 
agnetic systems results in characteristic temperatures lying in the range of 
0.01-0.1 K. For superparamagnetic nanoparticles (for which care has been 
taken to avoid direct contact between the particles by, e.g., applying a surfac- 
tant to a ferrofluid), exchange interaction and other interaction mechanisms 
can usually be discarded so that the dipolar interaction is the only relevant 
interparticle interaction. In addition, the size of the typical magnetic mo- 
ment (S ~ 10 2 — 10 5 magnetic spins) shifts the relevant temperatures up to 
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Figure 1: Magnetic energy versus 9 in the case of a longitudinal field for 
different values of the reduced field h = H/Hk- 



the range of a few kelvins, making it possible to observe effects of dipolar 
interaction in conventional magnetization experiments. 

The dipolar field, created by all other spins, at the position fi of the spin 
Si, is given by 

# = £?E<v«i-. ( 12 ) 

3 

where the term j = i is omitted from the summation, a is defined in such a 
way that a 3 is the mean volume around each spin, and 

Gij = -y (Sfijfij - 1) , (13) 

r ij 

r ij = r i ~ r j-i ?ij = — , (14) 

where 1 is the unit tensor. 

By introducing the dimensionless coupling constant 

_ (15) 

Cd - 4vra3 k B T ' (15j 

and noting that the dipolar energy = ^ J2i^j 'in • Hi, we can write the 
total dimensionless Hamiltonian of an interacting nanoparticle system as 

-(in = <7 £(3 • Uif + £ -h)+^Yl LOij , (16) 

i i i>j 
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where uj-ij = Sj • Gjj • Sj. Note that the interaction strength can also be 
measured by the temperature independent coupling parameter 

/ld= 2^ = 4^ C ' (17) 

which is the magnitude of the field, measured in units of the anisotropy 
field Hk-, produced at a given position by a dipole at a distance a. Here, 
c = V/a 3 is the volume concentration of particles. 

Dipole-dipole interaction is long-ranged and anisotropic, which makes 
it cumbersome to treat both analytically and numerically. Randomness in 
particle positions and anisotropy directions yields frustration and magnetic 
disorder leading to glassy dynamics for strongly interacting systems [18]. We 
will therefore use analytical and numerical treatment only for weakly inter- 
acting particle systems, while strongly interacting systems will be discussed 
in terms of collective spin glass behavior in section 3. 



2.2 Thermal Equilibrium Properties 

The thermal equilibrium average of any quantity B(s\, . . . , sjy) is given by 

{B) = z~J dTB ex p(-^) ' ( 18 ) 

where Z = JdT exp(— (37i) is the partition function and dT = YlidXli, with 
the solid angle dQi = d 2 Si/2ir. In the case of noninteracting spins, Eq. (18) 
has already been solved analytically for different quantities and anisotropies 
(for a review see, e.g., Ref. [8]). For isotropic spins, the magnetization is 
given by the Langevin function and the linear susceptibility follows a 1/T de- 
pendence, Xiso = ^P^om 2 , while the nonlinear susceptibility follows a 1/T 3 

dependence, xj^ = — g^Vo 7 ™ 4 - Magnetic anisotropy generally causes devi- 
ations from these well-known laws [8, 19-21] and also dipole-dipole interpar- 
ticle interaction will cause deviations from the Langevin behavior [22,23]. 
It is important to know the nature of such deviations in order to avoid 
confusion with, for example deviations due to quantum effects [24-26] . 

In the case of dipole-dipole interaction, the calculation of any thermody- 
namic property becomes a many-body problem and approximations and/or 
numerical simulations are needed in order to solve Eq. (18). Here we will 
use thermodynamic perturbation theory to expand the Boltzmann distribu- 
tion in the dipolar coupling parameter — an approximation that is valid for 
weakly coupled spins. 

The analysis presented in this section is restricted to spins with Hamil- 
tonians having inversion symmetry \H{fh) = H(—m)}. This assumption is 
valid for any kind of anisotropy or dipolar interaction, but if a bias field is 
applied in addition to the probing field, the condition of inversion symmetry 
breaks down. 
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2.2.1 Thermodynamic Perturbation Theory for Weakly Interact- 
ing Superparamagnets 

We will consider dipolar interaction in zero field so that the total Hamilto- 
nian is given by the sum of the anisotropy and dipolar energies TC = E a + Ed- 
By restricting the calculation of thermal equilibrium properties to the case 
£d <C 1, we can use thermodynamical perturbation theory [27, 28] to expand 
the Boltzmann distribution in powers of £d- This leads to an expression of 
the form [23] 

W = W & (l + S d F 1 + ±&F 2 + ---) , (19) 

where F\ is linear in and F 2 is (up to) quadratic in E^, while 

W a = Z' 1 exp(-/3£ a ) , (20) 

is the Boltzmann distribution of the noninteracting ensemble. Expressions 
for F\ and F 2 in are given in Appendix A.l. By keeping all averages weighted 
with W a , the thermal-equilibrium quantities calculated with this method will 
be exact in the anisotropy and only perturbational in the dipolar interaction. 
An ordinary high-temperature expansion corresponds to expanding Eq. (20) 
further in powers of (3. 

All results obtained below with the thermodynamic perturbation the- 
ory are limited to the case of axially symmetric anisotropy potentials (see 
Appendix A. 2), and all explicit calculations are done assuming uniaxial 
anisotropy (see Appendix B). 

2.2.2 Linear Susceptibility 

The equilibrium linear susceptibility is, in the absence of an external bias 
field, given by 

X=l^tf>. * = E(*-*)> (21) 

where h is a unit vector along the probing field direction and s z is the field 
projection of the net magnetic moment. Calculating (si) using thermody- 
namic perturbation theory yields an expansion of the susceptibility of the 
form 

X = ^ (ao + £dai + ffioa) , (22) 

with the general expressions for the coefficients a n given in Appendix A. 3. 
Simplified expressions for the coefficients can be obtained for some orienta- 
tional distributions of the anisotropy axes, such as parallel anisotropy axes 
and randomly distributed axes. 

For systems with parallel axes (e.g., single crystals of magnetic molec- 
ular clusters or a ferrofluid frozen in a strong field), the coefficients for the 
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longitudinal response read 



1 + 2S 2 

oo.ll = 3 ( 23 ) 

l + 452 + 45| 

1,11 ~ 9 ' ^ ' 

Kll = ~ 1 + A % + [(1 ~ S 2 ) (K-S) + 3S 2 (T - U)] (25) 
7 + 105 2 - 35Sf + 185 4 



+ 



315 



(1 - S 2 ) (K-K) + 3S 2 (T - \K) 



where C, 1Z (TZ), S, T and IA are lattice sums whose properties are discussed 
in Sec. 2.2.5. The properties of Si(a) are discussed in Appendix B. 

A common experimental situation is an ensemble of nanoparticles with 
the anisotropy axes oriented randomly (e.g., frozen ferrofluids). To obtain 
the susceptibility when the anisotropy axes are distributed at random, we 
average the general expressions for the a n over n, getting 

(26) 
(27) 

- Si) (K-K). (28) 
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Of), rand 


3 






G+rand 


5« 


1 _ _ 
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2 u 2,rand — 


~27 



1 

45 



Note that in the limit of isotropic spins (where Si — > 0) , the results for coher- 
ent axes and for random anisotropy duly coincide and agree with ordinary 
high-temperature expansions. 

It can easily be shown that ao, ra nd is independent of anisotropy for any 
type of anisotropy, not only axially symmetric [21]. Changing the coordinate 
system of Eq. (21) to the local one determined by the anisotropy direction of 
each spin, the field becomes a randomly distributed vector and by performing 
random averaging with respect to h by means of Eq. (84), one obtains 

x = J^m\^- Sl) / = 3k^T> (29) 

which is same expression as for isotropic spins. It is, however, only for the 
linear susceptibility term that randomly distributed anisotropy axes erase 
all traces of the anistropy. For the nonlinear susceptibility, the anisotropy 
is of importance even for systems with randomly distributed anistropy axis 
[8,19-21]. 
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2.2.3 Specific Heat 



The specific heat at constant volume can be obtained directly from the 
partition function 

g =/3 *Jl(l n *, = ^(,„*), (30) 

where, to take the a derivative, the coupling parameter ^ is expressed as 
£d = 2ahd [Eq. (17)]. As in the calculation of %, we consider only the 
zero-field specific heat. In that case, the term linear in ^ vanishes and the 
expansion of the specific heat to second order in ^ reads 

' = <A) + ^6 2 , (31) 



Nk B 

where the zeroth-order coefficient 

&o = ^(185 4 -355 2 2 + 10,S 2 + 7) , (32) 
315 

gives the specific heat in the absence of interaction [8]. The general formula 
for b 2 is given in Appendix A. 4. Again, it is possible to obtain simplified 
formulae for coherent anisotropy axes and for random anisotropy. In the 
first case (n, = n, Vi), we obtain 

62,11 = \{l-Sl-AaS 2 S' 2 -<j 2 [S 2 S'i + (S' 2 f}}n 
+\{2S 2 {l-S 2 ) + AaS' 2 {l-2S 2 ) 

+ a 2 {S'i-2[S 2 S l i + (S' 2 ) 2 ]}){1Z-TZ) 
+{s 2 + AaS 2 S' 2 + a 2 [S 2 S' 2 ' + (S' 2 ) 2 }}t , (33) 

where S' 2 = dS 2 /da. For randomly distributed axes, on averaging the general 
expression for b 2 over ft, one simply gets 

6 2 ,ra„d = \n. (34) 

This is the same correction term as that obtained for isotropic spins by 
Waller [29] and van Vleck [30] using ordinary high-temperature expansions. 



2.2.4 Dipolar Fields 

We are interested in calculating thermodynamical averages of the dipolar 
field, to introduce them in the expression for the relaxation rate in a weak but 
arbitrary oriented field, in order to obtain an expression for the relaxation 
rate of weakly interacting dipoles (we will argue in Sec. 2.3.3 that the effect 
of weak dipolar interaction can be accounted for by the thermodynamic 
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averages of the dipolar field). Because of the inversion symmetry of the 
anisotropy, only the square of the field will enter the low-field expression 
for the relaxation rate and not the field itself. In addition, the effects of 
longitudinal and transversal fields will be different. To be able to calculate 
the dipolar field is also of interest in the study of quantum tunneling of the 
magnetization of molecular magnets (e.g., Fes an d Mni2) [3,31] 

The general expressions for (£?m) and = — to second 

order in £d are given in Appendix A. 5. If we consider infinite systems, the 
index i can be removed since all spins have the same surroundings. For a 
system with aligned anisotropy axes, the averaged fields are given in terms 
of the lattice sums by 

<£f> = f [(1-S 2 )H + 3S 2 T\, 

<£l> = *f[3K-(l-S 2 )K + 3S 2 (K-K-T)] , (35) 
while for randomly distributed anisotropy axes they read 

(f) = f^ M> = §2^. (36) 

Note that the result for random anisotropy is identical to the result for 
isotropic dipoles. 



2.2.5 The Lattice Sums 

An essential element of the expressions derived for c v , and the dipolar 
fields are the following "lattice sums:" 1 

C = ^EE^ G ^ (37) 

i j+i 

*> = |EE^ 6 (38) 
K = ^EE^'-G^A (39) 

i j^i 

S = ^EEE^G^-G^ (40) 

r = ^EE(^-g^) 2 (41) 

i j^i 

U = ^EEE(^ G r#-G 3 v^) (42) 



1 h should be replaced by ft in the formulae for c„ and the dipolar fields. 



14 



The long range of the dipole-dipole interaction leads to a sample shape 2 
dependence of the physical quantities dependent on an external magnetic 
field [32,33]. In the expressions obtained for the susceptibility, this sam- 
ple shape dependence is borne by the slowly convergent lattice sums C, S, 
and U. If we consider "sufficiently isotropic" lattices, in the sense of fulfill- 
ing J2( r x) n = 2~X r 2/)™ = XX r z)™j such as cubic and completely disordered 
lattices (incidentally, the type of arrangements for which in the classical 
Lorentz cavity-field calculation the contribution of the dipoles inside the 
"small sphere" vanishes) these lattice sums vanish for large spherical sam- 
ples. The sums TZ, 7Z, and T, on the other hand, contain r^ 6 , which make 
them rapidly convergent and sample shape-independent. For sufficiently 
symmetric lattices 71 = 71 = 16.8, 14.5, 14.5 for simple cubic (sc), body- 
centered cubic (bcc), and face-centered cubic (fee) structures, respectively 
[34], and T = 13.54 (sc), 3.7 (bcc), 4.3 (fee). Note that since K - 71 = 0, 
some terms vanish in the above expressions for x, c v , and the dipolar fields. 

2.2.5.1 Sample Shape and Anisotropy Dependence The shape de- 
pendence of the linear susceptibility is illustrated in Fig. 2. The susceptibil- 
ity is calculated for small systems using both thermodynamic perturbation 
theory and a Monte Carlo technique [35], taking the dipolar interaction 
into account without any approximation. It can be seen that \ obtained 
by thermodynamic perturbation theory accurately describes the simulated 
susceptibility at high temperatures, while the results start to deviate at the 
lowest temperatures displayed. An estimate of the lowest temperature at- 
tainable by the thermodynamic perturbation theory is £d ~ 1/6, which is 
milder than the a priori restriction ^ <C 1. 

For the linear susceptibility, the zero-field specific heat as well as the 
dipolar fields, the anisotropy dependence cancels out in the case of randomly 
distributed anisotropy axes (at least for sufficiently symmetric lattices). In 
other cases the anisotropy is a very important parameter as shown for the 
linear susceptibility in Fig. 3 for an infinite (macroscopic) spherical sample. 
The susceptibility is divided by Xiso = ^Pl^om 2 in order to single out effects 
of anisotropy and dipolar interaction. It can be seen in this figure that the 
effect of the dipolar interaction is much stronger for a system with parallel 
anisotropy axesthan for a system with random anistropy (which coincide 
with the case of isotropic spins). At low temperatures, the susceptibility of 
a system with parallel anisotropy approaches that of Ising spins (calculated 
with an ordinary high-temperature expansion), while at high temperatures 
it is close to that of isotropic spins. 

The specific heat for uncoupled spins does not depend on the orientations 

2 Sample shape refers to the shape of the whole ensemble of nanoparticles not to the 
shape of the individual particles. The linear susceptibility exhibits a sample shape depen- 
dence, while the zero-field specific heat and the dipolar fields do not. 
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Figure 2: Equilibrium linear susceptibility in reduced units \ = x(HK/m) 
versus temperature for three different ellipsoidal systems with equation 
x /a 2 + y 2 /b 2 + z 2 /c 2 < 1, resulting in a system of N dipoles arranged 
on a simple cubic lattice. The points shown are the projection of the spins 
to the xz-plane. The probing field is applied along the anisotropy axes, 
which are parallel to the z axis. The thick lines indicate the equilibrium 
susceptibility of the corresponding noninteracting system (which does not 
depend on the shape of the system and is the same in the three panels); thin 
lines show the susceptibility including the corrections due to the dipolar in- 
teraction obtained by thermodynamic perturbation theory [Eq. (22)]; the 
symbols represent the susceptibility obtained with a Monte Carlo method. 
The dipolar interaction strength is h<i = £d/2c = 0.02. 
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Figure 3: Equilibrium linear susceptibility (x/Xiso) versus temperature for 
an infinite spherical sample on a simple cubic lattice. The dotted lines 
are the results for independent spins, while the solid lines show the results 
for parallel and random anisotropy calculated with thermodynamic pertur- 
bation theory, as well as for Ising spins calculated with an ordinary high- 
temperature expansions. We notice in this case that the linear susceptibility 
for systems with random anisotropy is the same as for isotropic spins calcu- 
lated with an ordinary high-temperature expansion. The dipolar interaction 
strength is h d = £ d /2cr = 0.004. 



of the anisotropy axes; however, the corrections due to the dipolar coupling 
do, as can be seen in Fig. 4). As for the linear susceptibility, the effect 
of dipolar interaction is stronger in the case of parallel anistropy than for 
random anistropy. 

2.3 Dynamic Properties 

At high temperatures, a nanoparticle is in a superparamagnetic state with 
thermal equilibrium properties as described in the previous section. At low 
temperatures, the magnetic moment is blocked in one potential well with 
a small probability to overcome the energy barrier, while at intermediate 
temperatures, where the relaxation time of a spin is comparable to the 
observation time, dynamical properties can be observed, including magnetic 
relaxation and a frequency-dependent ac susceptibility. 

For applications, such as magnetic recording media, it is necessary to 
know how different parameters will affect the relaxation time in order to 
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Figure 4: The specific heat per spin versus temperature for noninteracting 
spins (thick line) and weakly interacting spins with randomly distributed 
anisotropy axes (dashed lines) and parallel axes (thin lines) arranged on 
a simple cubic lattice. In each case, hd = £d/2cr = 0.003 and 0.006 from 
bottom to top. The inset shows the specific heat for noninteracting spins 
over a wider temperature interval. 



avoid spontaneous data erasure (caused by thermal fluctuations) on the 
lifetime of the device. Because of the ongoing effort to increase the infor- 
mation/volume ratio, it is of special importance to know how the dipolar 
interaction of densely packed nanoparticles will affect the relaxation time. 

In 1963 Brown [5] derived the Fokker-Planck equation for the probabil- 
ity distribution of the spin orientation, starting from the stochastic Gilbert 
equation, and calculated the relaxation time for particles with uniaxial 
anisotropy in a longitudinal field. More recent work on spins with nonaxially 
symmetric potential revealed a large dependence of the relaxation time on 
the damping coefficient A in the medium-weak damping regime [36-38] . Ex- 
periments on individual nanoparticles analyzed with accurate asymptotes of 
the relaxation time [39] gave damping coefficients in the range A « 0.05 — 0.5. 
Nonaxially symmetric potentials are, for example, created when applying a 
field at an oblique angle to a uniaxial spin. This oblique field can either be 
a bias field [40] or a nonlinear probing field [41]. In the case of interacting 
particles a transverse field component arises from the dipolar field of the 
surrounding particles. This explains the dependence on A of the blocking 
temperature that was first observed in numerical simulations by Berkov and 
Gorn [13] (see Fig. 5). 

The importance of including the damping in models describing the dy- 
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Figure 5: The out-of-phase component of the ac susceptibility versus temper- 
ature for two different values of the damping A = 0.1 and 1.0. ^M s uj = 0.03, 
^dem = M^V/k B , and h d = (/x /4tt)(c//3) with = 2.0. From Ref. [13]. 
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namic response of spins with nonaxially symmetric potentials (e.g., interact- 
ing uniaxial spins) tells us that models based only on how the energy barriers 
change [42, 43] necessarily overlook the precession of the magnetic moment 
(A — > oo) and therefore cannot account for the numerical results of Berkov 
and Gorn. A simple approach to include the damping in the dynamics of 
weakly coupled spin systems was proposed in Ref. [44]. 



2.3.1 The Equation of Motion 

We will present the equation of motion for a classical spin (the magnetic 
moment of a ferromagnetic single-domain particle) in the context of the 
theory of stochastic processes. The basic Langevin equation is the stochastic 
Landau-Lifshitz(-Gilbert) equation [5,45]. More details on this subject and 
various techniques to solve this equation can be found in the reviews by 
Coffey et al. [46] and Garcfa-Palacios [8]. 



2.3.1.1 Deterministic equations The motion of a magnetic moment 
can be described by the Gilbert equation [47] 

1 dm A _ dm 

— tt = mxH cfi m x — 43 

7 at 7m dt 

where 7 is the gyromagnetic ratio (which includes /j,q), A is a dimensionless 
damping coefficient, and the effective field is given by 

H eS = (44) 
on /dm 

The first term on the right side of Eq. (43) represents the precession of the 
magnetic moment about the axis of the effective field, while the second one 
is the damping term, which rotates rh toward the potential minima and is 
responsible for the dissipation of the energy. 

The Gilbert equation can be cast into the Landau-Lifshitz form [48] 

1 - d § = rn,H eS -^rnx(rn,H eS ) (45) 
7 at m 

with a "renormalized" gyromagnetic ratio 7 — ► 7/(1 + A 2 ) [8,46]. In the 
case of uniaxial anisotropy and a Hamiltonian given by Eq. (7), H e & = 
(H x /m)(fn ■ n)n + H , where Hk is the anisotropy field and H is an external 
field. The ferromagntic resonance frequency to for the precession about H c r 
is given by to = 7/i #cfr [49]. 



2.3.1.2 Stochastic Equations At T / the magnetic moment will 
interact with the microscopic degrees of freedom (phonons, conducting elec- 
trons, nuclear spins, etc.). The complexity of this interaction allows an 
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idealization, namely, to introduce them through a stochastic model. The 
simplest model is the Brownian, in which the interaction of m with the 
surroundings is represented by a randomly fluctuating magnetic field. This 
fluctuating field is necessarily combined with a dissipation (damping) term, 
and these two terms are linked by fluctuation-dissipation relations [50]. 

In the work of Brown [5] and Kubo and Hashitsume [45] the starting 
equation is the Gilbert equation (43), in which the effective field is increased 
by a fluctuating field yielding the stochastic Gilbert equation. This equation 
can, as in the deterministic case, be cast into the Landau-Lifshitz form as 

= ™ x [Has + h(t)\ -—mx{mx [H eS + b a (t)}} , (46) 

known as the stochastic Landau-Lifshitz-Gilbert (LLG) equation. The fluc- 
tuating field is assumed to be Gaussian distributed white noise 

<6fl,a(*)> = 0, (b Q , a (t)b^(t')) = 2D5 al3 5(t - t') , (47) 

with a, f3 = x, y, z. Garcia-Palacios and Lazaro [51] showed that the stochas- 
tic Landau-Lifshitz-Gilbert [Eq. (46)] and the simpler stochastic Landau- 
Lifshitz (LL) equation, 

--7r = mx[H cS + b {i (t)} m x (m x H cS ) , (48) 

7 at m 

both give rise to the same Fokker-Planck equation, describing the average 
properties of the magnetic moment, but with different Einstein- type rela- 
tions between the amplitude of the fluctuating field and the temperature: 

D LLG = -^ k -^, D LL = X k -^. (49) 
1 + A z 7m jm 

2.3.2 Relaxation Time in a Weak but Arbitrary Field 

We are interested in knowing how the relaxation time of uniaxial spins is 
affected by a weak field at an arbitrary direction, since it will allow us 
to study how the superparamagnetic blocking is affected by a field. This 
field dependence of the relaxation time can be obtained by expanding the 
relaxation rate T = 1/r in powers of the field components. As the spins have 
inversion symmetry in the absence of a field, T should be an even function 
of the field components, and to third order it is given by 

r~r (i + C ||£f + cj_£i) , (50) 

where £y and respectively are the longitudinal and transversal compo- 
nents of the field [given in temperature units, see Eq. (8)] with respect to 
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the anisotropy axis. Tq is the zero-field relaxation rate, which for low tem- 
peratures (a > 1) is given by Brown's result [5] 

To = --^'V* , (51) 
td V 71 " 

where td = m/ {2^\k-QT) is the relaxation time of isotropic spins. The 
coefficient cm can be obtained by expanding the expression for V in the 
presence of a longitudinal field [5, 52] 

r«||,a = 0) = ^-^[(l + h)e-^ 2 + (l-h)e-^- h ^} 

* r (i + ±£f). (52) 

There is no general expression for the relaxation rate in the presence of a 
nonzero transversal field valid for all values of the relevant parameters [53] , 
but Garanin et al. have derived a low-temperature formula valid for weak 
transversal fields [38], which can be used to determine the coefficient c±: 

r(£||=0,£ ± )~r [l + |F(a)£i], (53) 
F(a) = 1 + 2(2a 2 e) 1 /(^ 2 ) 7 (i + _L , _L ) . (54) 

Here a = Act 1 / 2 and 7(0, z) = Jq dtt a ~ l e~ l is the incomplete gamma func- 
tion. It can be noted that for the axially symmetric potential with a lon- 
gitudinal field, the only dependence on A is the trivial one in td, while in 
the nonaxially symmetric potential obtained with a transversal field the re- 
laxation rate will strongly depend on A through F(a), which is plotted in 
Fig. 6. 

Collecting these results, we can finally write the expression for the re- 
laxation rate in a weak field [44]: 



rv 



2.3.3 Relaxation Time of Weakly Interacting Nanoparticles 

The relaxation time for weakly interacting nanoparticles with uniaxial anisotropy 
can be obtained by inserting the thermodynamical averages of the dipolar 
fields (calculated in section 2.2.4) in the expression for the relaxation rate 
in a weak field [Eq. (55)], yielding T ~ T [l + \ (^j) + \F{a) (£[)]. Other 
models [42, 43, 54] are energy-barrier-based and therefore neglect the A de- 
pendence on relaxation time. For instance, the model by M0rup and Tronc 
[43] is basically the same as the one presented here in the particular case 
of high damping and random anisotropy. The M0rup-Tronc model predicts 
a decrease of the blocking temperature with increasing interaction strength 
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Figure 6: F versus a = \a x l 2 (solid line) and the overdamped value F = 1 
(dashed line). 



for weak interaction as was observed in high-frequency Mossbauer experi- 
ments, while the Dormann-Bessias-Fiorani model [42] predicts an increase 
of the blocking temperature with increasing interaction strength as com- 
monly observed in magnetization measurements. These discrepancies led to 
some controversy [55,56]. Berkov and Gorn [13] showed, by numerical inte- 
gration of the stochastic LLG equation [Eq. (46)], that for strong anisotropy 
(or weak interaction) the blocking temperature decreases with interaction 
(see Fig. 5), while for weak anisotropy (or moderate-high interaction) the 
energy barriers are governed by the interaction and hence grow with hd- 
An increase of the apparent blocking temperature is clearly the case for 
the strongly interacting nanoparticle samples (see Section 3) in which the 
relaxation time increases with due to spin-spin correlations (see Fig. 21). 

In order to determine the characteristics of the superparamagnetic block- 
ing we use the equilibrium susceptibility Xeq calculated using the thermody- 
namic perturbation theory Eq. (22) and the relaxation rate T obtained when 
the dipolar fields [Eq. (35) or (36)] are introduced in Eq. (55). Combining 
these expressions in a Debye-type formula 



X — Xeq 



T + ILO 



(56) 



provides us with a simple model for the dynamic response. The dynamic 
susceptibility of a large spherical sample with parallel anisotropy axis and 
a simple cubic structure is shown in Fig. 7. In the overdamped case, the 
blocking temperature is not noticeably affected by the dipolar interaction 
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while for low damping the blocking temperature decreases significantly as 
the interaction strength increases. These results are in agreement with the 
simulations by Berkov and Gorn [13] shown in Fig. 5. 
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Figure 7: Imaginary component of the dynamical susceptibility versus tem- 
perature (the real component is shown in the inset) for a spherical sample 
and spins placed in a simple cubic lattice. The anisotropy axes are all par- 
allel, and the response is probed along their common direction. The dipolar 
interaction strength h A = £ d /2a is h d = (thick lines), 0.004, 0.008, 0.012, 
and 0.016 from (a) right to left and (b) top to bottom. The frequency is 
lot d /(J = 2vr x 0.003. 

An interpretation of the strong damping dependence found in the pres- 
ence of a transverse field component was given in Ref. [41] — the transverse 
field creates a saddle point in the potential barrier. A thermally excited 
spin with high damping will "fall" directly back to the bottom of the poten- 
tial well if the thermal excitation is not sufficient for an overbarrier jump. 
On the other hand, a weakly damped spin in the same situation will precess 
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(~ 1/A times) about the anisotropy axis and therefore has an increased prob- 
ability for overbarrier jumps each time it passes close to the saddle point. 
In the case of noninteracting particles the transverse field component must 
come from either a nonlinear probing field [41] or a bias field [40], while the 
transverse field here naturally arises from the dipolar interaction. 

2.4 Numerical Methods 

Because of the long-range and reduced symmetry of the dipole-dipole inter- 
action analytical methods such as the thermodynamic perturbation theory 
presented in Sec. 2.2.1 will be applicable only for weak interaction. Nu- 
merical simulation techniques are therefore indispensable for the study of 
interacting nanoparticle systems, beyond the weak coupling regime. 

The Monte Carlo (MC) method can be used to efficiently calculate ther- 
mal equilibrium properties (c.f. Fig. (2). However, since it is an energy- 
barrier-based method, it will fail to generate dynamic features such as the 
precession of the spins, and will be able to generate the dynamic magneti- 
zation in the overdamped limit (A — > oo) only if an appropriate algorithm 
is used [35]. 

Using a Langevin dynamics approach, the stochastic LLG equation [Eq. (46)] 
can be integrated numerically, in the context of the Stratonovich stochas- 
tic calculus, by choosing an appropriate numerical integration scheme [51]. 
This method was first applied to the dynamics of noninteracting particles 
[51] and later also to interacting particle systems [13] (c.f. Fig. (5)). 

Because of the long-range nature of the dipolar interaction, care must 
be taken in the evaluation of the dipolar field. For finite systems the sums 
in Eq. (12) are performed over all particles in the system. For systems with 
periodic boundary conditions the Ewald method [57-59], can be used to 
correctly calculate the conditionally convergent sum involved. However, in 
most work [12, 13] the simpler Lorentz-cavity method is used instead. 

3 Strongly Interacting Nanoparticle Systems — 
Spin-Glass-Like Behavior 

A large number of applications use densely packed magnetic nanoparticles. 
It is thus important to know how interparticle interaction affects the physical 
properties of magnetic nanoparticle systems. In particular, it is important 
to understand how the thermal stability of magnetic recording media is 
changed by interparticle interactions due to the current effort of shrinking 
the volume of storage devices. 

It has been suggested that dense nanoparticle samples may exhibit glassy 
dynamics due to dipolar interparticle interaction [18]; disorder and frustra- 
tion are induced by the randomness in the particle positions and anisotropy 
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axis orientations. In order to investigate such systems, one needs to use the 
experimental techniques (protocols) developed in studies of spin glasses. Ex- 
amination of the effects of dipolar interactions using standard experimental 
protocols (zero-field-cooled ZFC / field-cooled FC magnetization) indicates 
no dramatic change in these quantities, and one can be misled to believe 
that the only effect of the dipolar interaction is to increase the blocking 
temperature due to enhanced energy barriers. However, glassy dynamics 
has been observed in strongly interacting nanoparticle systems [60-62] and 
for strongly interacting systems with narrow size distributions evidence has 
been given for a spin-glass- like phase transition [63-66]. 

In frozen ferrofluids the strength of the dipolar interaction between the 
single-domain nanoparticles can be continuously varied by changing the par- 
ticle concentration. With increasing particle concentration the magnetic be- 
havior may evolve from superparamagnetic to spin-glass-like. We will begin 
this section by recalling some fundamental properties of spin glasses. Fur- 
thermore, experimental results will be presented on a ferrofluid of Fe\- x C x 
nanoparticles. The glassy dynamics of dense nanoparticle samples will be 
compared with those of an Ising and a Heisenberg spin glass. 



3.1 Spin Glasses 

This section is intended as a brief introduction to spin glasses, focusing on 
the most recent work. For reviews on spin glasses, see, for example, Refs. 
[67-69]. 



3.1.1 Material 

The canonical spin glass consists of a noble metal (Au, Ag, Cu, or Pt) diluted 
with a transition metal ion, such as Fe or Mn. The magnetic interaction in 
such systems is mediated by the conduction electrons, leading to an indi- 
rect exchange interaction — the RKKY (Ruderman and Kittel [70], Kasuya 
[71], and Yosida [72]) interaction, whose coupling constant J(R) oscillates 
strongly with distance r between the spins as 

T f \ t cos(2fc F r + <gp) 

J(r) = Jo ^3 • ( 57 ) 

Here Jo and ipo are constants and kp is the Fermi wavevector of the host 
metal. Since the spins are randomly placed in the host metal, some spin- 
spin interaction will be positive and favor parallel alignment while other will 
be negative, thus favoring antiparallel alignment. 

The pure RKKY interaction is isotropic, and the canonical spin glass 
systems are therefore often referred to as Heisenberg spin glasses. However, 
some anisotropy is also present in those systems originating from dipolar 
interaction and interaction of the Dzyaloshinsky-Moriya (DM) type [73]. 
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The latter is due to spin-orbit scattering of the conduction electrons by 
non-magnetic impurities and reads 

Ebm = -Dij ■ (si x Sj), Dij ocfj x fj , (58) 

where is a random vector due to the randomness of the spin positions fj. 
The dipolar interaction is, as discussed in Secion 2.1.4, weak for atomic spin 
systems, while the DM interaction is enhanced by the presence of nonmag- 
netic transition-metal impurities [73,74]. However, for macroscopic spins 
(magnetic moments), the dipolar interaction is important and if the par- 
ticles are dispersed in a nonconductive medium (e.g., a frozen ferrofluid), 
it is the dominating interparticle interaction. If the magnetic particles are 
randomly placed, the dipolar interaction will be both positive and negative. 

A widely studied model system for an Ising spin glass is single crystals 
of Fe^Mni^TiOs with x « 0.5 [75-77]. Both FeTi0 3 and MnTi0 3 are 
antiferromagnets having the easy axis along the hexagonal c axis of the 
ilmenite structure. The Fe 2+ spins in FeTiC>3 are coupled ferromagnetically 
within a c layer and antiferromagnetically between adjacent c layers. In 
MnTiC>3, on the other hand, both the intralayer and interlayer coupling of 
Mn 2+ spins are antiferromagnetic. The compound Fe x Mni_ x Ti03 behaves 
as an Ising spin glass for 0.4 < x< 0.57, due to the mixing of ferromagnetic 
and antiferromagnetic interaction [77]. 

There are essential differences between the systems we have presented 
here, even in the limit of very strong anisotropy. The RKKY spin glasses 
are Heisenberg systems with random unidirectional anisotropy. Ferroflu- 
ids frozen under zero field are Heisenberg systems with random uniaxial 
anisotropy, while an Ising system is characterized by parallel uniaxial anisotropy. 

3.1.2 Spin Glass Models 

The Hamiltonian of an Ising spin glass, given by Edwards and Anderson 
(EA) [78], is 

N 

H = ~\^JijSiSj - H^Si, (59) 

i,j i=l 

where the spin Sj = ±1 and the coupling constants Jij are chosen from some 
random distribution fulfilling J2i,j Jij = in the case of a symmetric spin 
glass. In the EA model the spin-spin interaction is only of the nearest- 
neighbor type. The Sherrington-Kirkpatrick (SK) model [79] is the infinite- 
range version of the EA model. It is most useful as a basis for mean-field 
calculations. One such solution is the replica symmetry breaking theory of 
Parisi [80-82]. 

For a ferromagnet the order parameter is the magnetization, while for 
antiferromagnets it is the sublattice magnetization. For spin glasses the 
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magnetization is zero at all temperatures, and an appropriate order param- 
eter was proposed by Edwards and Anderson [78] as the average value of 
the autocorrelation function 

q E A = lim (Si(0) ■ . (60) 

The order parameter susceptibility, which diverges at the transition temper- 
ature, is the nonlinear susceptibility X2, defined as M/h = xo + X2h 2 + • • ■, 
where xo is the linear susceptibility [83]. The divergency of the nonlinear 
susceptibility was first shown on a Au(Fe) spin glass in 1977 by Chikazawa 
et al. [84] and more recently for a strongly interacting nanoparticle system 
by Jonsson et al. [64]. 

3.1.3 Critical Dynamics 

Close to the transition temperature T g , the dynamics of a spin glass system 
will be governed by critical fluctuations, but critical fluctuations are also of 
importance on experimental timescales quite far from T g . At temperatures 
both below and above T g , length scales shorter than the coherence length of 
the critical fluctuations 

£~L \e\- u , (61) 

will be dominated by critical fluctuations. Here Lq is a microscopic length 
scale and the reduced temperature e = (1— T/T g ). The coherence length can 
be transformed into a timescale according to conventional critical slowing 
down; the critical correlation time is given by 

t c ~ r m (£(T)/L ) 2 ~ T m \e\- zu , (62) 

where r m is a microscopic timescale. For spin glasses, r m ~ 10~ 13 s is 
the fluctuation time of an atomic moment. For nanoparticles, r m can be 
assigned to the superparamagnetic relaxation time of a single particle of 
average size, which, in the relevant temperature range for our studies, can 
be approximated by the Arrhenius-Neel expression [Eq. (6)]. 

At temperatures below T g , there is a crossover between critical dynamics 
on short length (time) scales and activated dynamics on long length (time) 
scales. The length scale of critical dynamics as a function of temperature 
is illustrated in Fig. 8. For T > T g , the system is in equilibrium on length 
(time) scales longer than £ (r c ) and hence the magnetic response is param- 
agnetic. The "freezing" temperature (Tf(r c )) of the crossover from a para- 
magnetic response to slow spin glass dynamics, as a function of observation 
time, can be obtained from experiments. Such a dynamic scaling analysis 
is performed in Sec. 3.2.2 for two samples of a nanoparticle system with 
two different volume concentrations. Dynamic scaling analyzes have given 
evidence for critical slowing down in a wide range of spin glass materials 
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Figure 8: Illustration of length (time) scales in spin glasses. 



[67,68,85]. The critical exponents obtained are rather scattered. However, 
the critical exponents of Ising and Heisenberg systems are clearly different 



3.1.4 Nonequilibrium Dynamics 

Aging phenomena in glassy materials were first discovered and thoroughly 
investigated in the field of structural glasses [88]. Magnetic aging in spin 
glasses was first observed by Lundgren et al. in 1983 [89]. It was found that 
the ZFC relaxation depends on the wait time t w during which the system 
has been allowed to age at the measurement temperature, before applying 
the probing field and recording the magnetization as a function of time 
t. The measurement protocol is illustrated in Fig. (9) and ZFC relaxation 
measurements on a Ag(ll at% Mn) spin glass are shown in Fig. 10. It can be 
seen that the relaxation rate S(t) = h~ 1 dM/dlogt exhibits a maximum at 
t ~ t w [89, 90]. Many different types of materials were later shown to exhibit 
aging and nonequilibrium dynamics, including polymers [91], orientational 
glasses [92], gels [93], and ceramic superconductors [94]. 

Spin glass systems represent ideal model systems for studying nonequilib- 
rium dynamics experimentally, numerically, and theoretically. The dynamic 
properties of spin glasses can accurately be investigated by superconductive 
quantum interference device (SQUID) magnetometry on specially designed 
magnetometers with low-background field (< 1 mOe) and optimized tem- 



[86,87]. 
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Figure 9: Schematic representation of a ZFC relaxation experiment. The 
sample is cooled to the measurement temperature T m under zero field, after 
a wait time t w , a small probing field h is applied and the ZFC magnetization 
is recorded as a function of time. 

perature control (AT < 100 fiK) [95]. Spin glasses exhibit nonequilibrium 
dynamics on all timescales from r m to infinity. Experimental studies are, 
however, limited to a finite time window tjr m ~ 10 8 — 10 17 , while numerical 
simulations are limited to much shorter timescales t/r m ~ 1 — 10 5 . Because 
of the larger spin-flip time of a magnetic moment compared to an atomic 
spin, interacting nanoparticle systems can be investigated in an experimental 
time window in between those of ordinary spin glasses and numerical sim- 
ulations. In addition, the strength of the dipolar interaction can be tuned 
by the particle concentration of the nanoparticle system. Strongly interact- 
ing nanoparticle systems are therefore also interesting as model systems for 
glassy dynamics. 

3.1.4.1 The Droplet Model The "droplet" theory [96-100] is a real- 
space theory, based on renormalization group arguments for the Ising EA 
model with a continuous distribution of independent exchange. It makes 
predictions concerning the nonequilibrium dynamics within the spin glass 
phase. Important concepts are domain growth — growth of the coherence 
length for equilibrium spin-glass order, and temperature chaos. Both these 
concepts can also be applied to less anisotropic spin glasses [87, 101-103] 
and generalized to other glassy systems [104]. 

At each temperature the equilibrium spin glass state is considered to 
consist of a ground state plus thermally activated droplet excitations of 
various sizes. A droplet is a low-energy cluster of spins with a volume L d 
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Figure 10: (a) ZFC susceptibility versus time on a logarithmic scale for the 
Ag(ll at% Mn) spin glass. T m = 30 K and h = 0.5 G. (b) Relaxation rate 
S(t) = h~ l J^- t of the ZFC susceptibility curves shown in (a). 



and a fractal surface area L ds . The typical droplet free-energy scales as 

/ r \0 

F ty P _ T(T) ^_ j ^ T(T) - Je 9v , (63) 

where 6 is the stiffness exponent and T(T) is the stiffness. The droplet free 
energy is broadly distributed, and because of the presence of configurations 
that are almost degenerate with the ground state, the distribution of Fl will 
have weight down to zero energy: 

PL(F L ) w P{Fl 'S\ p(0)>0. (64) 

If 9 < 0, large droplets can be flipped at a low energy cost so that the large 
droplets will not be stable against small fluctuations and the system will be 



31 



paramagnetic. Hence, a negative value of 9 indicates that the system is below 
its lower critical dimension [96, 105]. On the other hand, if 9 > 0, very few of 
the large-scale droplets will be thermally activated since F L yp > k^T. Since 
Pl(Fl) has non-zero weight near zero energy, a certain fraction of droplets 
will be thermally active and dominate most of the equilibrium physics. 

The dynamics of droplets is considered to be a thermally activated pro- 
cess. The energy barrier for annihilation of a droplet will scale as 

4 yp ~A(T)(A^, A(T)~Je^, (65) 

where A(T) sets the free-energy scale of the barriers and ip is an exponent 
satisfying 9 < ip < d — 1. The characteristic time tl that a thermally 
activated droplet will last for is given by an Arrhenius law 
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Lh (66) 



where tq{T) is the unit timescale for the activated process. For activated 
hopping processes the unit timescale is not simply given by the real micro- 
scopic timescale [106], which is r m ~ %/J ~ 10~ 13 s in spin systems. A 
plausible choice for tq(T) is instead the critical correlation time r c [Eq. (62)] 
as proposed in Ref. [100]. The Arrhenius law implies that droplets of length 
scale L = Lx(t) 



L T (t) 



k B Tln(t/T (T)) 



(67) 



A(T) 

can be activated within a timescale t. 

Lritw) will be the characteristic length scale of equilibrium spin glass 
order at a time i w after a quench from a temperature above T g to a temper- 
ature T in the spin glass phase. 

In a magnetization measurement the system is probed by applying a 
small magnetic field. The magnetization arises through the polarization of 
droplets. Since this polarization also is a thermally activated process, it will 
affect droplets of size L(t), where t is the time elapsed after the application 
of the magnetic field in a ZFC-relaxation experiment or t = 1/uj in an 
ac experiment at a given angular frequency to. The peak observed in the 
relaxation rate of the ZFC magnetization at t ~ t w (see Fig. 10) has, within 
the droplet model, been interpreted as a crossover from quasiequilibrium 
dynamics at L(t) < L(t w ) to nonequilibrium dynamics at L(t) > L(t w ) (see 
Ref. [107] for a detailed discussion). 

According to the droplet theory, typical spin configurations of a pair 
of equilibrium states at two different temperatures, say, T\ and T2, are es- 
sentially the same on short length scales much below the so-called overlap 
length Lat, but completely different on large length scales much beyond 
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Lat- This temperature chaos is due to a subtle competition between en- 
ergy and entropy in the spin glass phase. In the limit of small temperature 
differences \AT/J\ <C 1, the overlap length between the two temperatures 
Ti and T 2 = T x + AT is supposed to scale as [99, 100, 108] 



AT 



J 



C 



(4 - 20) 



(68) 



where £ is the chaos exponent. 

Experimentally, temperature chaos can be evidenced by aging the system 
at a temperature Tj and changing the temperature to a T 8 + AT. If the 
temperature shift is large enough, the equilibrium domain configurations at 
Ti and Tj + AT are completely different on the length scales relevant for the 
experiment. Hence, if the magnetic response at Tj + AT is the same as after 
a direct quench, the system appears rejuvenated. 



3.1.4.2 Experiments: Aging, Memory and Rejuvenation After 
discovery of the aging effect in the spin glass phase, experimental protocols 
with temperature steps and cyclings were proposed [109]. These experiments 
showed not only rejuvenation effects but also that spin glass order, character- 
istic of different temperatures, can coexist on different length scales; hence 
the spin glass keeps a memory of its thermal history. 

T + 

T 2 




' 1 . J I . j 

t w i t w 2 t 

Figure 11: The experimental procedure of a "memory" experiment with two 
halts on cooling. 

A simple experimental protocol [85] was employed to illustrate memory 
and rejuvenation effects [110-113]; the sample is cooled from a high temper- 
ature with one (or more) halts of the cooling at one (or more) temperatures 
in the spin glass phase. The experimental procedure of a double-stop exper- 
iment is illustrated in Fig. 11. The ac susceptibility is subsequently recorded 
on heating. An "ac memory" experiment is shown in Fig. 12 for a Ag(ll 
at% Mn) spin glass and in Fig. 13 for an Feo.5oMno.5oTi03 spin glass. The 
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ac susceptibility is measured on cooling with two intermittent stops. Dur- 
ing a stop the ac susceptibility relaxes downward as shown in the inset of 
Fig. 12. The level of the ac susceptibility is hence related to the age of the 
system (a lower susceptibility indicates an older system). As the cooling 
is resumed, the ac susceptibility merges (quite rapidly) with the reference 
curve, and the system is rejuvenated. On subsequent reheating the ac sus- 
ceptibility shows a dip around each aging temperature — the system has 
kept a memory of each equilibration at constant temperature. The memory 
experiment is an efficient tool to study spin-glass-like properties in various 
materials [114-118]. 




T(K) 

Figure 12: x"( w ) versus temperature for the Ag(ll at% Mn) sample mea- 
sured during cooling (circles) and during the subsequent reheating (pluses). 
Two intermittent halts were made during the cooling: one at 27 K for 10, 000 
s and another at 23 K for 30, 000 s. The susceptibility measured on con- 
stant cooling and on constant heating is shown as reference (solid lines). 
The arrows indicate the cooling and heating curves, respectively. The inset 
shows \"{oj) versus time during the halt at T = 23 K; uj/2tt = 510 mHz. 
The cooling and heating rate is ~ 0.2 K/min. 

For the Ag(ll at% Mn) sample, the ac susceptibility curve measured on 
cooling lies below the curve subsequently measured on heating, except close 
to the lowest temperature (see Fig. 12). If the aging at different tempera- 
tures is accumulative, the heating curve would appear older than the cooling 
curve and therefore lower in amplitude. For the Feo.5oMno.soTi03 sample 
this is indeed the case (see Fig. 13). The nonaccumulative behavior observed 
in the cooling-heating curves of the ac susceptibility of the Ag(ll at% Mn) 



34 




T(K) 

Figure 13: a) x" versus temperature measured on cooling (diamonds), with 
two intermittent stops at 20 K for 10, 000 s and at 16 K for 30, 000 s, and 
on the subsequent reheating (pluses), for the Feo.5oMno.5oTi03 sample. The 
reference cooling and heating curves are drawn with solid lines. Inset: x" ~ 
Xr e f versus temperature derived from the heating curves for single and double 
stops. w/2-7r = 510 mHz. 

sample can qualitatively be explained by rejuvenation during cooling and 
heating due to strong temperature chaos. For the Feo.5oMno.5oTi03 sample, 
the double-memory experiment shown in Fig. 13 indicates that temperature 
chaos does exist, but that it is much weaker than for Ag(ll at% Mn) . 

The memory experiments presented above yield, in a simple and illus- 
trative way, information about aging, memory, and rejuvenation effects for 
the two spin glass systems. However, cooling and heating rate effects as 
well as memory and rejuvenation phenomena are all mixed in a nontrivial 
way. Specially designed thermal protocols have been used to quantitatively 
investigate aging and rejuvenation proprieties of the Feo.5oMno.5oTi03 and 
Ag(ll at% Mn) sample [87, 102, 103, 119]. It was shown that the tempera- 
ture dependence of Lr(t) is stronger for the Ag(ll at% Mn) sample than 
for the Feo.5oMno.5oTi03 sample. Also, the rejuvenation effects are much 
stronger for the Ag(ll at% Mn) sample. 

3.2 FeC Nanoparticle Systems 

We now focus on a ferrofluid of single-domain particles of the amorphous 
alloy Fe\- x C x (x « 0.2-0.3). The particles were coated with a surfactant 
(oleic acid) and dispersed in a carrier liquid (xylene). The particle shape is 
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nearly spherical (see Fig 14) and the average particle diameter d = 5.3 ± 
0.3 nm. The saturation magnetization was estimated to M s = 1 x 10 6 A/m, 
the microscopic flip time to To = 1 x 10~ 12 s, and the uniaxial anisotropy 
constant K = 0.9 x 10 5 J/m 3 [65]. The interparticle interaction can be 
varied by changing the particle concentration of the ferrofluid. The strength 
of the interaction for a given concentration is determined by the anisotropy 
constant and the saturation magnetization according to Eq. (17) with the 
parameters given above hd « 0.56c, where c is the volume concentration of 
nanoparticles. The samples studied here originate from the same batch as 
those in Refs. [65, 114, 120, 121]. Earlier studies use samples from different 
batches having slightly different physical properties [63, 64, 122-124]. 

Figure 15 shows the real and imaginary part of ac susceptibility versus 
temperature for three different particle concentrations of the FeC sample: 
c = 0.06, 5, and 17 vol%. With increasing concentration, the peak in the 
ac susceptibility is shifted to higher temperatures and the curve is simul- 
taneously suppressed. This behavior is different from that of the weakly 
interacting nanoparticle systems shown in Figs. 5 and 7. In this section we 
will argue that the dynamics of the 5 and 17 vol% samples is spin-glassdikc, 
and hence fundamentally different from the superparamagnetic behavior of 
noninteracting and weakly interacting nanoparticle systems. 




Figure 14: TEM picture of typical Fei- x C x nanoparticles. 
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Figure 15: Ac susceptibility vs temperature at frequencies lo/2it = 125 Hz 
(filled symbols) and u>/2ir = 1000 Hz (open symbols). 

3.2.1 Nonequilibrium Dynamics 

Magnetic aging can be evidenced by measuring the ZFC relaxation at con- 
stant temperature after a fast cooling through the transition temperature 
using different wait times before applying the magnetic field. The experi- 
mental procedure was depicted in Fig. 9 together with a typical measurement 
on a Ag(ll at% Mn) spin glass in Fig. 10. In a noninteracting nanoparticle 
system, the low-field ZFC relaxation is governed only by the distribution of 
relaxation times of the particles and their temperature dependence. It does 
not depend on the wait time at T m before applying the probing field, 3 as 

3 The only wait time dependence that could exist is the adjustment of the position of 
the magnetic moment to the Boltzmann distribution at T m . This adjustment does however 
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Figure 16: S'(i) versus time on a logarithmic scale for the 5 vol% sample 
obtained from ZFC relaxation measurements with i w = 300 s (open symbols) 
and 3000 s (filled symbols); h = 0.05 Oe. 



was shown experimentally in Ref. [60]. The relaxation rate S(t) of the ZFC 
magnetization measured at different temperatures between 20 and 40 K are 
shown in Fig. 16 for the 5 vol% sample. The measurements are repeated 
for two different wait times 300 and 3000 s. A clear difference between the 
S(t, t w ) curves for t w = 300 and 3000 s can be seen for all temperatures 
< 40 K, presenting evidence for glassy nonequilibrium dynamics at those 
temperatures. The shapes of the S(T) curves are, however, rather different 
from those of canonical spin glasses (see Fig. 10). A peak in S(t) at t ~ t w as 
in ordinary spin glasses is observed only at T = 30 K. The difference arises 
largely from the strong temperature dependence of the individual particle 
relaxation time compared to the almost temperature independent relaxation 
time of individual spins. 

The 5 vol% sample has also been investigated with ac susceptibility mea- 
surements. x'(T) is shown in Fig. 17 for a low frequency. The insets show 
how the ac susceptibility relaxes with time if the cooling is halted at T s = 33 
or 23 K. The relative relaxation of the ac susceptibility is smaller than for 
ordinary spin glasses. %' relaxes more in absolute units than \" and x' nas 
therefore been chosen to illustrate the nonequilibrium dynamics. For a wait- 
ing time of 30, 000 s at low temperature Ax" /Xref ~ 3 % compared to « 8 



% for the Feo.5oMno.5oTi03 sample (see Fig. 13) and 
at% Mn) sample (see Fig. 12). 



30 % for the Ag(ll 



only involve intrawell rotation and is therefore a much faster process than the slow cooling 
to T m and the time needed to stabilize the temperature. 
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Figure 17: x' versus temperature for the 5 vol% sample. The insets show 
how x' relaxes with time if the cooling is halted at 33 or 23 K. / = 510 mHz 



Memory experiments with one and two temporary stops during cooling 
are shown in Fig. 18. The nonequilibrium effects are more clearly revealed 
by subtracting the reference curves obtained on constant cooling and reheat- 
ing. The features are qualitatively similar to those of ordinary spin glasses 
(see Sec. 3.1.4.2); during a halt in the cooling, the system ages — the ac 
susceptibility decreases. When the cooling is resumed, the ac susceptibility 
slowly regains the reference level. On the subsequent heating, the ac sus- 
ceptibility exhibits a "dip" centered around T s in a single-stop experiment. 
In a double-stop experiment it exhibits two dips if the two halts are well 
separated in temperature [Fig. 18(a)], but if the two halts are close in tem- 
perature [Fig. 18(b)], only one large dip is observed on heating as a result of 
the two aging processes. In both cases, the difference curve x' ~ x' TC { °f the 
double-stop experiment equals the sum of x' ~ x' rc f °f the two single-stop 
experiments. As for the Feo.soMno.soTiOs sample, the ac susceptibility mea- 
sured on heating lies below the one measured on cooling. In conclusion, the 
rejuvenation effect in particle systems is weaker than in ordinary spin glasses 
and as a consequence, the memory of both one and two aging processes is 
better preserved on reheating than in the spin glass case. 

This nanoparticle sample exhibits strong anisotropy, due to the uniaxial 
anisotropy of the individual particles and the anisotropic dipolar interaction. 
The relative timescales (t/r m ) of the experiments on nanoparticle systems 
are shorter than for conventional spin glasses, due to the larger microscopic 
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Figure 18: x'(T) — Xref(^) versus ^ measured on cooling (symbols) and 
heating (lines) for the 5 vol% sample. Circles — the cooling was halted at 
33 K for 5400 s; diamonds — the cooling was halted at (a) 23 K for 36 000 s, 
(b) 28 K for 25, 200 s; pluses— the cooling was halted at T x = 33 K for 
t Wl = 5400 s and at (a) T 2 = 23 K for t W2 = 36 000 s, (b) T 2 = 28 K for 
t W2 = 25 200 s. Xrcf is shown in Fig. 17. lj/2tt = 510 mHz. 

flip time. The nonequilibrium phenomena observed here are indeed rather 
similar to those observed in numerical simulations on the Ising EA model 
[125, 126], which are made on much shorter time (length) scales than exper- 
iments on ordinary spin glasses [127]. 

3.2.2 A Spin Glass Phase Transition? 

Aging and nonequilibrium dynamics indicate but give by no means evi- 
dence for a thermodynamic phase transition at finite temperature to a low- 
temperature spin-glass-like phase. For example, two-dimensional (2D) spin 
glasses (T g = 0) have been shown to exhibit similar nonequilibrium dynam- 
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ics as 3D spin glasses [128, 129]. A detailed analysis of the magnetic response 
close to the assumed transition temperature is needed in order to evidence a 
spin glass phase transition (see, e.g., Ref. [85] and references cited therein). 

The existence of a second-order phase transition can be evidenced from 
critical slowing down [Eq. (62)] approaching the phase transition from the 
paramagnetic phase. Defining a criterion to determine the freezing temper- 
ature Tf associated with a certain relaxation time, it is possible to derive 
the "transition" line (r c ) between thermodynamic equilibrium and critical 
dynamics as in Fig. 8. In ac susceptibility measurements the relaxation time 
r = I/to, and a possible criterion for the freezing is where x"(T, o>) attains 
a certain fraction, say, 15%, of its maximum value. 




T(K) 

Figure 19: ac susceptibility versus temperature for the 17 vol% sample. 
w/2vr = 0.017, 0.051, 0.17, 0.51, 1.7, 5.1, 17, 55, 170 Hz. 
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Ac susceptibility data for a large set of frequencies are shown in Fig. 19 
for the 17 vol% sample and in Fig. 20 for the 5 vol% sample. These curves 
were used to extract the T{(l/u}) plotted in Fig. 21. The corresponding 
blocking temperatures for the 0.06 vol% sample are plotted in the same fig- 
ure as a reference for the behavior of a noninteracting system of the same 
particle ensemble. A dynamic scaling analysis according to critical slow- 
ing down [Eq. (62)] with the microscopic timescale given by an Arrhenius 
law [Eq. (6)], r c (T f ) = r exp(^/fc B 7f)(T f /T g - l) zv , was performed for the 
two concentrated samples. Two assumptions concerning the anisotropy en- 
ergy was used: (i) ^4 = 0, which correspond to a temperature-independent 
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Figure 21: Relaxation time r = uj~ 1 versus Tf . For the 5 vol% and 17 vol% 
samples the lines are fits to the critical slowing down relation [Eq. (62)] with 
the parameters given in Table 1. The assumptions E = and E = 500 yield 
exactly the same line. For the 0.06 vol% sample Tf is the superparamagnetic 
blocking temperature defined as the maximum of x" • 



microscopic flip time; and (ii) A/ks = 500 K, which is approximately the 
anisotropy barrier energy for a particle of average size. The values obtained 
for zu, T g , and ro in each case are given in Table 1. The quality of the fits 
of the experimental data to Eq. (62) are equally good for assumptions (i) 
and (ii). In fact, in Fig. 21, the line corresponds to both assumptions; in the 
experimental temperature-frequency interval the two assumptions cannot 
be distinguished. In addition, the values of T g and the critical exponents 
depend quite strongly on the criterion used when determining Tf . The error 
bars on the exponents obtained are therefore large. 

In a full dynamics scaling analysis, the imaginary component of the 
dynamic susceptibility of a spin glass is scaled according to [130] 



X"(?>) _„ 



Xeq(T) 



e"G(uJT c ), T > T g (69) 



where u = 1/t and G(x) is a scaling function. The asymptotic behavior of 
G{x) oc x y with y = 1 and (3/ zu for small and large values of x, respectively. 
It is shown in Fig. 22 that the x" data for the 17 vol% sample could be 
collapsed into a master curve according to this relation. In this figure the 
assumption A = is used, but an equally good collapse could be obtained 
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Table 1: Parameters obtained from dynamic scaling analysis 



Sample (vol%) 


E/k B (K) 


zv 


T g (K) 


T o (s) 


17 





11.4 


48.8 


2 • 10" 8 


17 


500 


8.8 


49.9 


5 • lO" 11 


5 





10.3 


36.0 


2 • 10~ 5 


5 


500 


6.4 


37.9 


1 • io- 8 



by assuming A = 500 K (and changing e@ to r m ). The critical exponent 
(3 = 1.0 ± 0.3. For the 5 vol% sample, the dynamic susceptibility could not 
be scaled according to Eq. (69). In addition, if data for larger values of 
r (obtained from ZFC relaxation data) are included in the critical slowing 
down analysis, deviations from the power law is observed [65]. To conclude, 
the dynamic scaling analysis indicates that the 5 vol% sample does not 
exhibit a thermodynamic phase transition although it clearly exhibits spin 
glass dynamics. Because of the temperature dependence of r m , a static 
scaling analysis as performed in Ref. [64] is a crucial additional tool to 
disclose a possible spin glass phase transition in interacting nanoparticle 
systems. 

The time dependence of the dynamic correlation function q(t) was inves- 
tigated numerically on the Ising EA model by Ogielski [131]. An empirical 
formula for the decay of q{t) was proposed as a combination of a power law 
at short times and a stretched exponential at long times 

q(t) = cT x e- wtv , (70) 

where c, x, w, and y are temperature-dependent parameters. q{t) follows 
a pure power-law behavior below T g . It was shown in Ref. [132] that 
[Xeq — x'(t)]/Xeq measured on the Feo.5oMno.5oTi03 spin glass sample also 
behaves according to Eq. 70. In Fig. 23 [x C q — x'(*)]/Xeq is shown for the 
17 vol% sample at temperatures around T g . The temperature dependence 
of the exponent x is shown in the inset. The behavior of the 17 vol% sam- 
ple shown here is similar to that observed in the numerical simulation on 
the EA Ising model [131] and to the experimental result obtained for the 
Feo.5oMno.5oTi03 spin glass [132]. However, the stretched-exponential be- 
havior is less pronounced than for the Feo.5oMno.5oTi03 spin glass. One 
reason for this is that the investigated timescales t/r m are shorter for the 
nanoparticle sample because of the larger value of r m . In the numerical 
simulation, performed on even shorter timescales, it was possible to observe 
the stretched-exponential behavior by investigating q(t) far from T g . The 
resolution of our experimental data does not allow such an investigation. 



44 




Figure 22: Scaling of e'^x" '(T,u) /xeq data for T > T g for the 17 vol% 
sample. The assumption A = is used and (3 = 1.0. The other parameters 
are those of Table 1. The two lines are the asymptotic behavior G(x) oc x 
for small values of x and G(x) oc x^l zv for large x. Inset: critical slowing 
down analysis on a log-log scale. 

3.2.3 Dynamics in a Field 

All measurements presented so far are performed with a probing field in 
the linear response regime and without an external bias field. We will now 
investigate the effect of a nonzero bias field on the nonequilibrium dynamics 
(still probing the system with a weak magnetic field in the linear response 
regime) . The question as to whether a spin-glass-like phase exists under a fi- 
nite field in a strongly interacting nanoparticle system will not be addressed. 
Figure 24 shows the ac susceptibility as a function of temperature with bias 
dc fields in the range - 250 G for 1 and a 5 vol% samples. The magnetic 
field affects the superparamagnetic behavior of individual spins (as discussed 
in section 2) as well as the nonequilibrium dynamics. It can be seen in the 
figure that the dynamic response of the two systems in low fields is rather 
different, while with increasing bias field the dynamic response of the two 
systems becomes more and more similar. This indicates that the effects of 
dipolar interaction are suppressed by a sufficiently strong magnetic field. 

The effect of an external field on the glassy dynamics can be studied 
by recording the ac susceptibility as a function of time for different bias 
fields, x" normalized by its value just after the quench, x"(^o)> is shown in 
Fig. 25 for the 5 vol% sample. It can be seen in this figure that the ac relax- 
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Figure 23: [x C q — x'(t)]/Xeq vs time for the 17 vol% sample. Inset: temper- 
ature dependence of the exponent x in Eq. (70). 



ation diminishes with increasing bias field. For the highest field, almost no 
relaxation exists. The field has in addition the effect that it makes the mea- 
surement more noisy. The effects of the field is also temperature-dependent 
and the ac susceptibility is more affected by a field at high temperatures. A 
qualitatively similar result was found for the (Feo.i5Nio.85)75Pi6BeAl3 spin 
glass [121]. 

3.3 Discussion 

We have seen that the magnetic properties of a strongly interacting nanopar- 
ticle system are of spin-glass-like nature, and hence very different from the 
superparamagnetic behavior of noninteracting systems and weakly interact- 
ing spins discussed in section 2. Any model for interparticle interaction 
based on a modified superparamagnetic behavior [42-44, 54] will therefore 
fail to describe the dynamics of a strongly interacting nanoparticle system. 

The value of the dipolar coupling parameter hd [Eq. (17)] determines 
the strength of the dipolar interaction, but the width of the distribution 
of energy barriers is equally important for the dynamic properties. For 
example, Jonsson et al. studied [133] a 7-Fe2C>3 nanoparticle sample that 
had a value of £d comparable to the 17 vol% sample investigated here, but 
with a much broader energy barrier distribution. That sample showed glassy 
dynamics, but it did not exhibit a spin glass phase transition. 

It would be interesting to examine the nonequilibrium dynamics of a 
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Figure 24: ac susceptibility versus temperature for different superimposed 
dc fields; H = 0, 50, 100, 150, 200, 250 Oe. u/2tt = 125 Hz. 

strongly interacting system for different concentrations (different dipolar 
coupling strengths). Such a study could give more detailed insight into how 
the aging and rejuvenation phenomena depend on the interaction strength 
and the experimental timescale. 

4 Summary and Conclusion 

The effects of interparticle dipolar interaction in magnetic nanoparticle sys- 
tems have been investigated. For weak dipolar coupling strengths equilib- 
rium quantities can be calculated using thermodynamic perturbation theory, 
treating the anisotropy exactly and the dipolar interaction perturbatively. 
Such an approach is always valid at sufficiently high temperatures, but it 
is valid around the blocking temperature only in the case of weak dipolar 
coupling. In a simple model, the relaxation rate is modified by the dipolar 
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Figure 25: ac susceptibility versus time for different superimposed dc fields; 
H = 0, 40, 80, 120, 160, 200 Oe. w/2vr = 510 Hz. 



interaction in the same way as by a field, where the field components of the 
dipolar field can be calculated using thermodynamic perturbation theory 
[44]. The dipolar field also plays an important role for quantum tunneling 
of the magnetization of crystals of magnetic molecules [3] . 

Superparamagnetic blocking depends strongly on the damping parame- 
ter in the case of weak-medium damping due to the transverse component of 
the dipolar field [38] . Any energy-barrier-based model overlooks that damp- 
ing dependence and can therefore only be valid in the overdamped case. The 
effect of a magnetic field is to decrease the relaxation time. Hence, in the 
case of weak interparticle interaction, a decrease of the blocking tempera- 
ture is predicted and also observed in high-frequency measurements, such 
as Mossbauer spectroscopy [43] and Langevin dynamics simulations [13]. 
However, for stronger interparticle interaction, the dipolar interaction does 
not modify but creates energy barriers, and hence the blocking temperature 
increases with increasing interaction strength. A blocking temperature that 
increases with the interaction strength is commonly observed in magnetiza- 
tion measurements. 

The existence of glassy dynamics in dense frozen ferrofluids, specifically 
strongly interacting nanoparticle systems with randomness in the particle 
positions and anisotropy axes, has been evidenced by experimental tech- 
niques developed in the study of conventional spin glasses. Hence, the dy- 
namics of such systems is radically different from simple superparamagnetic 
blocking. The nonequilibrium dynamics observed in strongly interacting 
nanoparticle systems exhibits qualitatively similar aging, memory, and re- 
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juvenation effects as ordinary spin glasses, but the aging and rejuvenation 
effects are weaker. The differences observed can be explained at least qual- 
itatively by the longer microscopic relaxation time of a magnetic moment 
compared to an atomic spin. Only strongly interacting nanoparticle systems 
with a narrow anisotropy barrier distribution have been shown to exhibit a 
phase transition to a low-temperature spin glass phase. 
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A Thermodynamic Perturbation Theory 

A.l Expansion of the Boltzmann Distribution in the Dipolar 
Coupling Parameter 

Thermodynamic perturbation theory is used to expand the Boltzmann dis- 
tribution in the dipolar interaction, keeping it exact in the magnetic anisotropy 
(see Sec. 2.2.1). A convenient way of performing the expansion in powers of 
is to introduce the Mayer functions fij defined by 1 + fij = exp^u;^), 
which permits us to write the exponential in the Boltzmann factor as 

expi-PH) = exp(-/3£ a ) JJ(1 + f tj ). (71) 

i>j 

Expanding the product to second order in the fij gives 

na + fij) = i + e d G! + \iiG 2 + o(d), (72) 

i>j 

where [134] 

Gi = £u^, (73) 

i>j 

G 2 = + ^2^2^ij^klQik:jiqil:jk, (74) 

i>3 i>j k>l 

and the symbol qik-.ji annihilates terms containing duplicate pairs: qik-.ji = 
\{2 - 5 ik - Sj^il + 5 ik ){l + 5ji). 
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To obtain the average of any quantity B, we introduce the expansion (72) 
in both the numerator and denominator of (B) = JdTB exp(— (3Ti)/ JdT exp(— PTC), 
and work out the expansion of the quotient, getting 

(B) * (B) a + ^ (B Gi) a + ±&[(B G 2 ) a - (B) a <G 2 ) a ]. (75) 

Here we have utilized the fact that the single-spin anisotropy has inversion 
symmetry [W a (— Sj) = W a (sj)] in the absence of a bias field and that (Gi) a = 
J2i>j u ij ) = since a dipole does not interact with itself. 
We have now obtained expressions for F\ and F 2 in Eq. (19): 

Fi = Gi , F 2 = G 2 - (G 2 ) a . (76) 

To complete the calculation, we need to obtain averages of low-grade powers 
of s weighted by the noninteracting distribution (moments), which is the 
only place where one needs to specify the form of E a . In the next section 
we will do that for systems with axially symmetric anisotropy. 

A. 2 Averages Weighted with an Axially Symmetric Boltz- 
mann Factor 

Assuming an axially symmetric potential, the anisotropy energy of E a (s- n) 
will be an even function of the longitudinal component of the magnetic 
moment s ■ n. The averages we need to calculate are all products of the 
form I m = (Il™=i (c*n • s)) a , where the c n are arbitrary constant vectors. 
Introducing the polar and azimuthal angles of the spin (p), we can write 
I m as 

j = Jo 27 " d(f Jo (I'd sin g Un=i (4 • s) exp [-(3E a (s ■ n)] 
Jq W dip Jq d"9 sin "9 exp[— f3E a (s • ft)] 

For odd m, I m is an integral of an odd function over a symmetric interval 
and hence I m = 0. To calculate the susceptibility and specific heat to second 
order in £d> we require I 2 and I4, which will be calculated using symmetry 
arguments similar to those employed to derive the a = unweighted averages 
(see, e.g., Ref. [135]). 

Note that I 2 is a scalar bilinear in c\ and c 2 . The most general scalar 
with this property that can be constructed with the vectors of the problem 
(c*i, c*2, and ft) has the form 

I 2 = A ci ■ c 2 + B (ci ■ n)(c 2 • n) . 

To find the coefficients A and B, one chooses particular values for the c n : 

(i) If ci || c 2 _L n then I 2 = A. Thus, setting n = z and c\ = c 2 = x, one 
has s-n = cos ■d = z and (c\ -s)(c 2 -s) = (1 — z 2 ) cos 2 ip, so the integral 
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reads 

J 27r dip cos 2 <pj\ dz{\ - z 2 ) exp[-f3E a (z)} 



A 



J 27r d<pf\dzexp[-f3E a (z)] 
= *[1-<AJ = 



l-5 2 



where S 2 = (P 2 {z)) & is the average of the second Legendre polynomial 
-^2(2) = \ (3z 2 — l) over the noninteracting distribution. 

(ii) If ci || <?2 || n, then I 2 = A + B. Putting n = c\ = c 2 = z the integral 
is given by 

A + B= I-idzz 2 exp[-(3E a (z)] = 2 _ 1 + 2S 2 
[ ' . d~, exp| JI--A-) 



Therefore, since I 2 = ((ci ■ s){c 2 - s)) a , we get the following for the second- 
order moment: 

1 - S 2 

((ci • s)(c 2 • s)) a = — - — ci • c 2 + 5 2 (ci • n)(c 2 • n) . (77) 

We can similarly calculate I4 by constructing the most general scalar fulfill- 
ing certain properties, getting 

((ci • s){c 2 ■ s)(&3 ■ s)(c4 • s)) a = A 4 [(ci • c 2 )(c 3 • c 4 ) + (ci • c 3 )(c 2 • c 4 ) 

+(ci • c 4 )(c 2 • c 3 )] 

+A 2 [(ci • c 2 )(c 3 • n)(c 4 • ft) + (ci • c 3 )(c 2 • n)(c 4 • n) 
+(ci • c 4 )(c 2 • n)(c 3 • ft) + (c 2 • c 3 )(ci • n)(c 4 • n) 
+ (c 2 • c 4 )(ci • n)(c 3 • ft) + (c 3 • c 4 )(ci • ra)(c 2 • ra)] 
4-^(01 • n)(c 2 ■ n)(c 3 ■ n)(c 4 • ft), (78) 

where A 2 and A4 are combinations of the first Si(cr) 

A 2 = i(5 2 -5 4 ), A 4 = |-f + l. (79) 

Therefore, Eq. (78) involves S 2 as well as £4 = (P 4 (z)) a , the average of the 
fourth Legendre polynomial -P 4 (z) = | (35 z 4 — 30 z 2 +3) with respect to 

Finally, introducing the following tensor and scalar shorthands 

T = i^l + 5 2 nn, (80) 
A, _ „ „ A 2 



A = ^A 4 l + - 7 =p=nn, ft = S4 - 3— ^ , (81) 
' A 4 A 4 
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where 1 is the identity tensor, the results for the moments can compactly 
be written as 

((c 1 -s)(c 2 -s)) a = (c!-r.c 2 ) (82) 
((ci ■ s)(c 2 ■ s)(c 3 ■ s)(c 4 • s)) a = (ci ■ A • c 2 )(c 3 ■ A ■ c 4 ) 

+ (ci ■ A ■ c 3 )(c 2 • A • c 4 ) 

+ (ci-A-C4)(c 2 -A-c 3 ) (83) 

+ Q(ci ■ H)(c 2 • n)(c 3 • n)(c 4 • n) , 

which facilitates the manipulation of the observables. 

The quantities Si are calculated in the case of uniaxial anisotropy in 
Appendix B. Note finally that in the isotropic limit (Si — > 0), Eqs. (77) and 
(78) reduce to the known moments for the isotropic distribution [134, 135] 

((ci • s)(c 2 ■ s))iso = I ci ■ c 2 , (84) 
((ci ■ s)(c 2 ■ s)(c 3 ■ s)(c 4 ■ s))iso = js[(ci ■ c 2 )(c 3 ■ c 4 ) 

+(ci • c 3 )(c 2 • c 4 ) 
+ (ci ■ c 4 )(c 2 -c 3 )] • (85) 

These expressions are formally identical to those for the average of a quan- 
tity involving the anisotropy axes Hi, when these are distributed at random 
Jt J2i f(™i) ~^ 1 = /■ For instance, for arbitrary n-independent 
vectors v\ and v 2 , we have 

— ^2(vi ■ Hi)(v 2 ■ Hi) — ► (vi ■ H)(v 2 ■ H) = \ vi ■ v 2 . (86) 



A. 3 General Formulae for the Coefficients of Susceptibility 

The general expression for the equilibrium linear susceptibility is given by 
Eq. (22) with the following expressions for the coefficients 

«0 = i-rY, h - T i- h ( 8? ) 



N 
1 

N 



fl i = ^EEMIVG^-IV)./* (88) 

i j^i 



«2 = 4^^ / '' !r '' G ' / ' rrG; '' r ' ) ' / ' 

+ f EEE^ ( r * • <V r i ■ G jk -r k ) ■ h 

i jjki k^j 



i j^i 1 'ij 
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+ 2(h ■ Ai ■ hi? 



+ Q(h ■ n,i) 2 (ni ■ hjY 



+ S 2 



-jvEEC 1 - 1 '^) 



(h- Ai - h)(fij ■ Gij ■ Ai ■ dj ■ Hj) 
+ 2(h • Ai • • Hj) 2 
+ £l(h ■ Hi) 2 (Hi ■ Gij ■ Hj) 
1-S 2 



+ S 2 (fij ■ G^ ■ Ti ■ G^ ■ Hj) 



(89) 



where Gij, fij and fij are defined in Eq. (13), and T, A, and in Eqs. (80) 
and (81) and also involve the Si(a). 

When calculating these coefficients, the same type of averages appear as 
in the isotropic case (see Refs. [30, 134] for details of the calculation) and 
with the same multiplicities. The only difference is the weight function and 
hence the formulas required to calculate those averages [Eqs. (77) and (78) 
instead of Eqs. (84) and (85)]. 

A. 4 General Formula for the Coefficient b 2 of Specific Heat 

In the general expression (31) for the specific heat the coefficient 60 is given 
by Eq. (32), while b 2 reads 

Nb 2 = l{ 2 {l-S 2 )-AaS' 2 -a 2 S'i}Y,Y,n" 

i j^i 

+ \{2S 2 (l- S 2 )+4oS' 2 (l-2S 2 ) 
+ a 2 [S' 2 \l-2S 2 )-2(S' 2 ) 2 )} 

*EE^ 2 ' 21 



[nj ■ Hi) + (fij ■ Hj) ] 



i j^i 



+{S 2 2 + 4aS 2 S' 2 + a 2 [S 2 S' 2 ' + (S' 2 ) 2 }} EEfa ' ' n i? ( 90 ) 

i j^i 

where /' = df /da. General formulae for in the case of uniaxial anisotropy 
are given in Appendix B. 

A. 5 Dipolar Field 

The dipolar field averages we want to calculate are (£, 2 ^) = • Hi) 2 } and 
(£?_!_) = — (Cj 2 ||)- The general expressions for these quantities are [44] 

= f E^ 1 -^) (Hi -Gij -dj -Hi) 
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+3S 2 {Hi ■ Gij ■ nj) 2 } , (91) 

(£,±) = f E^ + ^^^-G^-n,) 

j 

- (1 - S 2 ) (fit ■ ■ Gij ■ Hi) 

-3S 2 {Hi ■ Gij ■ ntf] . (92) 

Since we want to use the dipolar fields in order to examine the blocking 
behavior (see Sec. 2.3.3), one may wonder about the validity of these field 
averages below the superparamagnetic blocking, where the spins are not in 
complete equilibrium. However, since at those temperatures the spins are 
still in quasiequilibrium confined to one of the two wells, we can repeat the 
derivation of the algorithm (77) restricting the phase space for integration to 
one well. In this case, averages of the form {s ■ vi) a do not vanish, and should 
be considered together with ((s ■ v\)(s- v 2 )) a , which, being even in s, is not 
modified. The extra terms associated with (s- V\) a will, however, vanish if 
the overall state is demagnetized, and Eqs. (91) and (92) are recovered. 



B Si for Uniaxial Anisotropy 



The thermodynamical average Si(a) over the Legendre polynomials Pi occur 
in the expressions for the susceptibilities, the specific heat, and the dipolar 
fields in Sec. 2.2. For uniaxial anisotropy these averages read 



J-l 

In particular, Sq = 1 and S 2 = \ (3z 2 — 1) can be written 

S 2 = 



J.) 

2a J 



(93) 



(94) 



The one-spin partition function Z a = j\ dz exp(crz 2 ) can be written in 
terms of error functions of real and "imaginary" argument as 



Z a 



v / vr/Verfi( v / o : ), a > 
V / 7r7H" erf (v / H) 5 ° < 



(95) 



The less familiar erfi(x) is related to the Dawson integral D(x), so in the 
easy-axis case one can write Z a = (2e CT / 'y/a)D(y/a) and compute D(x) with 
the subroutine DAWSON of Ref. [136]. 

For I > 2, the Si can be computed using the following homogeneous 
three-term recurrence relation [137]: 



2a 



(2/-l)(2/ + 3) 



Si- 



2a 



21 + 1 



I - 1 
21 - 1 



Si. 



1 + 2 
21 + 3 



Si+ 2 



, (96) 
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The derivative of any 5/ can be computed by means of the differential re- 
currence relation [23] 

, dSi (1-1)1 21(1 + 1) 

1 da (2Z-1)(2Z + 1) '~ 2+ 3(2Z- l)(2/ + 3) ' 

(2Z + 1)(2Z + 3) + 3 v ' 

The approximate behavior of 52 and 54 for weak (\a\ <C 1) and strong 
(\a\ S> 1) anisotropy are 

^cr + ^a 2 + • • • |cj| < 1 
= { + . (98) 

S4(a) = { + + ■■■ cr»l . (99) 

3/i i 5 i 35 



g(l + i + J r ) + ... 
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